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scalar particles which are coupled to Abelian gauge fields. Scalars may inter- 
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Q>^ ' act with gravity in a non-conformally invariant way and thus their fluctua- 

CN ] tions can be amplified during infiation. These fiuctuations lead to creation of 
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electric currents and produce inhomogeneous distribution of charge density, 
resulting in the generation of cosmological magnetic fields. We address the 
question whether these fields can be coherent at large scales so that they may 
seed the galactic magnetic fields. Depending upon the mass of the charged 
scalar and upon various cosmological (critical fraction of energy density in 

X" 

^ ' matter, Hubble constant) and particle physics parameters we found that the 



magnetic fields generated in this way are much larger than vacuum fiuctua- 
tions. However, their amplitude on cosmological distances is found to be too 
small for seeding the galactic magnetic fields. 
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I. FORMULATION OF THE PROBLEM 



The idea that our galaxy could possess a magnetic field dates back to the (independent 
but simultaneous) works of Fermi (motivated by the origin of high energy cosmic rays) 
and Schwinger (motivated by the origin of galactic synchrotron emission). Since then, 
a lot of work has been done both from the experimental and theoretical side. Listening to 
observations |^,^ large scale cosmological magnetic fields can be estimated from Faraday 
rotation effects and Zeeman splitting of (hyperfine) spectral lines. Listening to theory the 
main puzzle is connected with the dynamical origin of large scale magnetic fields. 

It is usually assumed that the observed magnetic fields were (exponentially) amplified 
by galactic dynamo mechanism from pre-existing seed magnetic fields, coherent over scales 
of the order of 100 pc at the time of galaxy formation. The amplitude of the necessary seed 
fields is quite uncertain and depends on many details of the dynamo mechanism and on 
cosmological parameters Typical vaules of the seeds range from 10^^^ to 10^^^ Gauss 

at the decoupling epoch. 

In general, two types of ideas are considered for the explanation of the seed magnetic 
fields. The first one is related to astrophysics. For instance, a Biermann battery mechanism 
can be postulated at the level of protogalaxy [|| and this will possibly lead to a mean 
current able to provide a source term for the evolution equations of the magnetic fields in 
the galactic plasma. The observation of large scale magnetic fields associated with clusters 
of galaxies also suggested the appealing possibility of connecting (inter) galactic magnetic 
fields with active galactic nuclei JlOl (in this picture the dynamo mechanism is not essential 
since the rotation of a cluster is much smaller than the rotation of a galaxy). The second 
type of ideas relies on the the interplay between particle physics and cosmology at different 
moments in the life of the Universe. In this framework various models were discussed such 
as cosmological defects [|lT|, phase transitions [1^ electroweak anomaly |13|, inflation []6|,p!4 



string cosmology |jT5[, temporary electric charge non-conservation ||T6[ or breaking of the 
gauge invariance M. 



2 



The main problem of most particle physics mechanisms of the origin of seed fields is 
how to produce them coherently on cosmological scales. All causal proposals (related, for 
example, to the bubble collisions at the phase transitions) may produce sufficiently large 
magnetic fields only on sub-horizon scales. Infiation, quite effective in producing density 
fluctuations at super-horizon scales, fails to amplify directly the vacuum fluctuations of the 
electromagnetic fleld because of its conformally invariant coupling to gravity. 

An attractive and very economical idea on the possible primordial origin of the galactic 
magnetic flelds was suggested in |^. In short, it is based on the following observation. While 
the coupling of electromagnetic fleld to the metric and to the charged flelds is conformally 
invariant (this is not necessarily true in the models with dynamical dilaton the coupling 
of the charged scalar fleld to gravity is not. Thus, vacuum fluctuations of the charged scalar 
fleld can be amplifled during inflation at super-horizon scales, leading potentially to non- 
trivial correlations of the electric currents and charges over cosmological distances. The 
fluctuations of electric currents, in their turn, may induce magnetic flelds through Maxwell 
equations at the corresponding scales. The role of the charged scalar fleld may be played 
by the Higgs boson which couples to the hypercharge fleld above the electroweak phase 
transition; the generated hypercharged fleld is converted into ordinary magnetic fleld at the 
temperatures of the order of electroweak scale. 

No detailed computations were carried out in in order to substantiate this idea. The 
suggestion of was further developed quite recently in [jl^ for the standard electroweak 



theory with an optimistic conclusion that large scale magnetic flelds can be indeed generated. 
In [|18| a supersymmetric model was studied. 



The aim of the present paper is to re-analyze this proposal and compute the amplitude 
and the spectrum of seed magnetic flelds that may be generated because of ampliflcation 
of zero-point fluctuations of the scalar fleld during inflation. Our set-up resembles the one 



of Ref. [0. We suppose that an inflationary phase was followed by a radiation dominated 
phase and we compute the charged particle production associated with the change of the 
metric. This allows to deflne the spectrum and magnitude of the current fluctuations at the 



beginning of the radiation era. Taking the current distribution as an initial condition, we 
study the plasma-physics problem of the relaxation of such an initial condition. We compute 
finally the magnetic fields, which survived possibly until the time of galaxy formation. 

Depending upon the mass of the charged scalar and upon various cosmological (critical 
fraction of energy density in matter, Hubble constant) and particle physics parameters we 
found that the magnetic fields generated in this way are much larger than vacuum fluctua- 
tions, in agreement with qualitative conclusions of refs. . However, in contrast with 



17| , their amplitude on cosmological distances is found to be too small, by many orders of 
magnitude, in order to seed the galactic magnetic fields. We trace back this difference in the 
conclusions to the discrepancy in the results obtained for the Bogoliubov coefficients (ap- 
pearing in the problem of scalar particle production) and to the treatment of the relaxation 
of electric currents in conducting media during the radiation dominated epoch. 

The plan of our paper is the following. In Section II we will discuss the amplification 
of the charged scalar field in an expanding cosmological background. In Section HI we will 
study the connection between the amplification of the charged scalar and the production 
of charge and current density fluctuations. In Section IV we will develop a curved space 
description of the evolution of charge and current fluctuations within a kinetic (Landau- 
Vlasov) approach. We will apply our analysis to the physical case of an ultra-relativistic 
plasma prior to decoupling. Section V contains some phenomenological applications of our 
formalism. We will mainly be concerned with the generation of large scale magnetic fields in 
various models of cosmological evolution and with the possible occurrence of charge density 
and current fluctuations at large scales. Section VI contains our concluding remarks. 

II. AMPLIFICATION OF A COMPLEX SCALAR FIELD DURING INFLATION 



In this paper we will only consider scalar electrodynamics. Possible generalizations of 
our results to theories containing more scalar fields (for example, for electroweak theory or 
its extensions) are straightforward. Since fermions are conformally coupled to gravity, their 
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gravitational production is too small to generate any substantial seed magnetic fields [|T7 . 

Consider the action of a (massive) charged scalar field minimally coupled to the back- 
ground geometry and to the electromagnetic field (hypercharge field, if the standard model 
is assumed): 

S = J (fx^[{V^(t)rV^(l) - - -^T^pT'^^ (2.1) 

where T)^ = — ieA^, T^y = d[^Au], g^u is the four-dimensional metric and g its determi- 
nant. 

We will suppose that the background geometry is described, in conformal coordinates, 
by a Friedmann- Robert son- Walker (FRW) line element 

g,iu = a^Vtiu, ds^ = a?'{vi)[drf - dx^], (2.2) 

where rj^i, is the usual Minkowski metric. 

It is convenient to introduce rescaled fields $ = a0 and = aA^j,. Correspondingly, 
we denote by E and B the electric and magnetic fluctuations in curved space. They are 
related to the usual fiat space fields £ and S by a time-dependent rescaling involving the 
scale factor: 

E = a^S, B = a^B. (2.3) 

In terms of the rescaled fields the action is 

S = J d^xdr][r]^''D^^*Dy(!> + ~ ^^^^'>^*^ " (2-4) 

where the prime denotes differentiation with respect to the conformal time coordinate r) and 
Df^i = — ieA^, F^y = d[i^iAy]. This is simply the action of electrodynamics in flat space- 
time with a time dependent mass term for the scalar fleld. From this form of the action 
it is obvious that there is no direct ampliflcation of electromagnetic flelds during inflation. 
Moreover, for conformal coupling of the scalar fleld to gravity the term proportional to a" 
is absent and the scalar particle production is supressed by the charged boson mass. 



To compute the magnetic field fluctuations we wiU use a perturbative approach. Namely, 
we will firstly compute the scalar particle production omitting completely the coupling of 
the scalar field to the gauge field (i.e. neglecting the back reaction, as it can be checked to 
be, a posteriori, self-consistent). 

The classical evolution equations for $ in the case = are simply given by: 



a" 



$"_V2$ <l> + mV$ = 0. (2.5) 

a 

We will often use also decomposition of the complex scalar field in terms of two real fields, 
$ = ($1 +i$2)/y2. 

Once the background geometry is specified, the amplified inhomogeneities in the field 
$ can be computed. Suppose, as in that the history of the Universe consists of two 
different epochs. A primordial phase, whose background evolution is not exactly known, 
and a radiation dominated phase where the scale factor 0(77) evolves linearly in conformal 
time. A continuous (and differentiable) choice of scale factors is then represented by 

V 111 J 

+ + (2.6) 

where a is some effective exponent parametrizing the dynamics of the primordial phase of 
the Universe. Notice that if a = 1 we have that the primordial phase coincides with a de 
Sitter inflationary epoch. In practice we will consider a = 1 or slightly deviating from it. 

As a result of the change in the behaviour of the scale factor occurring at —rji the modes 
of $ will be parametrically amplifled. Deflning xi = kr]i and fi = mrji we can write the 
Bogoliubov coefficients for a = 1 (standard inflation) and for cosmologically interesting case 
xi <^ 1 (long ranged fluctuations) and /i <^ 1 (small scalar mass, potentially giving rise to 
large scalar fluctuations): 



1 t X 1 

ak = e'8^7r<! \ — - + ^ — - + 

2r(|)^i 2r(f)/ii 

3 1 



1 ^^(z-l)/iV4 



Xi +C(/i3), 



2r(f)^i v^r(i) 

/j, = e-fv^( ""'J '''!\ + \ 1 + ^^^'^1 v^U(9(/.t). (2.7) 



More general expressions can be found in Appendix A where we also give details about the 
calculation. The terms kept in the expansion for small Xi and fi is such that the unitarity 
condition for the Bogoluibov coefficients is satisfied: lafcp — = 1. This property of the 
truncation is quite important: when discussing charge density fiuctuations we will see that 
interesting cancellations between the leading terms arise. 

In the opposite limit, for k > l/r]i the mixing coefficients are exponentially suppressed 
as exp[— AA;?7i] |19|. The coefficient A depends upon the details of the transition between the 
inflationary and radiation dominated phases. The existence of an exponential suppression 
in the mixing coefficients is quite important from a general point of view: it ensures gentle 
ultra-violet properties for the physical quantities we ought to compute. 

We should mention that the leading contribution to the amplification coefficients has 
been computed in different contexts pO| . However, we are not only interested in the leading 
behaviour of the Bogoliubov coefficients but also in the corrections whose contribution is 
relevant when the leading contribution cancels, as in the case of charge density fluctuations 
(see the following Section). 



III. CHARGE AND CURRENT DENSITY FLUCTUATIONS 

The Bogoliubov coefficients obtained in the previous Section specify the probability of 
charged particle creation. The fluctuations in the scalar fleld induce also fluctuations in the 
electric current associated with the U{1) symmetry of our action. The fluctuations in the 
charge and current density act as a source term for the evolution of the fluctuations in the 
gauge flelds. 

We will ignore, for the moment, the effect of the plasma conductivity which is essential 
for the calculation of induced magnetic flelds (it will be discussed in the following section) 
and we will simply consider the structure of the current-current correlators. In the curved 
space, the conservation equation for the current is given by: 



where 

f ^ie{(t)*di'(t)-(t)d^(t)*). (3.2) 
It is convenient to introduce the rescaled current 

J' = V^g'^'j.. (3.3) 

which can be expressed as 



(3.4) 



in terms of conformal fields. Notice that in this last expressions the index appearing in the 
derivatives is raised (or lowered) using the Minkowski metric and not the curved metric. 

If we average on the vacuum we have clearly that (J^) = 0. However, the fluctuations 
of the same quantity for two space-time separated points are not zero. By defining the two- 
point functions of the field operators $i and $2 

g{x, y) = {^,{x)^M) = ($2(x)$2(2/)), (3.5) 
the two-point function of the charge and current density fiuctuations can be written as 

{Mx)My)) = 2e'[g{x,y)^^g{x,y) - Q^^i^^y)Q^,^i^^y)}^ (3-6) 

where x = {x, rj) and y = {y, r) {rj and r are two different conformal times). In the case of 
the vacuum (no amplification took place) the two-point functions are simply 

In the case of flat space-time the current density fluctuations can be expressed as 

iJM)) - / / f - *.)^e-<-.<-) (3.8) 

When the background passes from an inflationary phase to a radiation dominated phase we 
have instead that the two-point functions of the field operators can be written as 
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f d?k 




1 (27r)3 


Ittfcl 



|2„*/ 



+ ak(il9k{r])9k{r) + al(3kgl{v)9*k{r) 



e-'^"; (3.9) 



where Qkiv) are given by Eqs. ( |A.12D and where r = x — y . 



In order to define properly the charge and current density fiuctuations at different scales 
it is helpful to introduce an averaging procedure both over space and over time. In the case 
of the charge density fluctuations we will deflne the charge fluctuations inside a patch of 
volume Lr" and within a typical time T as 

Qlr=^ J d'x J d'y Jdr]J dr( Jo(x, r/) Jo(y, r))iyL,T(f , - eo)W^L,r(y, r - ^o), (3.10) 

where Wi^r are the smearing functions selecting the contribution of the correlator inside a 
given space-time region. We can choose the smearing functions to be, for instance, the so- 
called top hat function which has a sharp edge. It is deflned as WL^r{.x,rj) = 1 for |x| < L, 
\ri\ < T and it is equal to zero otherwise. This choice is, however, not so useful for our 
case. In fact, one can show that the Fourier transform of the top hat proflle goes to zero, 
for large k, as {kro)~^. Unfortunately this behaviour (ultimately related to the shrap edge 
of the proflle) could create spurious effects in the charge and current density fluctuations. 
A better choice is, for our purposes, the gaussian smearing function 

WL,r{x,v) = (3.11) 

This expression smears the high frequency modes more efficiently than the top hat function. 
By inserting Eqs. ( |3.9| ) into Eq. ( |3.10| ) we obtain, after regulatization over time, in the limit 
L <^ m^^ and for > Vi 

2 _ 2 - 2 

QIt - J J d^y J d^k J d^p\ap/3k - a^Ppl'^e-'^-^ e'^'^+P^^^-^^ . (3.12) 

A similar procedure can be carried on in the case of the current density fiuctuation. 

As we will see in the next Section, the relevant quantity to be computed in order to 
estimate the size of the gauge field fiuctuations is given by the trace of 
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q2 

((V X J)fc(V X J)i) = e,kaejibg^^,{Mx,7])Jj{y,T)). (3.13) 



Therefore, the regularized quantity we are interested in is 

(V X J)l^ = ^jd^xj £y Jdr^l rfr((V x J),(V x J)k)WL,T{^,r] - ^o)WL,T{y,r- 6)- 

(3.14) 



After regularization over time we find, using Eq. (|3.9|) into Eq. ( p.l4|) in the hmit Tn^Q/rji <^ 
1, that for L ^ m^^ and for > Vi our correlator becomes 

(V X J)l^ ^ I d'x I dPy I d'k I d'pV{k,p)e-i-^e'^^+^^<''y\ (3.15) 

where 

V{k,p) = [kY - {k ■ p}^]\ak + (3k\^\ap + I3p\\ (3.16) 



IV. MAGNETIC FIELDS FROM CHARGE AND CURRENT DENSITY 

FLUACTUATIONS 

Once we know the charge and current fluctuations at the beginning of the radiation 
dominated epoch we can compute the induced fluctuations in the electromagnetic fields from 
the Maxwell equations. The complete solution of this problem is hardly possible in a general 
case because of the complicated (and model dependent) dynamics of the infiaton decay and 
of the reheating of the Universe after inflation. We shall use an extremely simplified picture 
of this process. Namely, we will suppose that the plasma right after inflation is in thermal 
equilibrium with some temperature T, up to small (at super-horizon scales) perturbations 
of the distribution functions leading to charge and current fluctuations. Though this picture 
may be, in general, far from reality we expect it to reproduce the physics of the magnetic 
field generation quite accurately, because the rate of equilibration of the electromagnetic 
reaction is large compared with the rate of the Universe expansion after inflation and may 
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be bigger than the rate of the inflaton decay. In fact, our result depends essentially on 
the density of the charged particles and their typical momentum only. Thus they may be 
qualitatively applicable even to situations with large deviations from thermal equilibrium. 

Within this set-up the problem can be related to the relaxation of an initial perturbation 
in the distribution function. Consider an equilibrium homogeneous and isotropic conducting 
plasma, characterized by a distribution function fo{p) common for both positively and neg- 
atively charged ultrarelativistic particles (for example, electrons and positrons) . Suppose 
now that this plasma is slightly perturbed, so that the distribution functions are 

U{x,p,r]) = fo{p) + 6U{x,p,r]), f-{x,p,r]) = fo{p) + 6f-{x,p,r]), (4.1) 

where + refers to positrons and — to electrons, and p is the conformal momentum. The 
Vlasov equation defining the curved-space evolution of the perturbed distributions can be 
written as 

+ V ■ + e[E + V B) ■ = [——] , (4.2) 



dr] dx dp V drj ) coW 

^ + ^.^-e{E + vxB)-^=(^) (4 3) 

dr] dx dp \ dr] /cou' 

where the two terms appearing at the right hand side of each equation are the collision 
terms. This system of equation represents the curved space extension of the Vlasov-Landau 
approach to plasma fluctuations |2l|j2^. All particle number densities here are related to 



the comoving volume. 



Notice that, in general v = p/yfP + mio?. In the ultra-relativistic limit v = p/\p\ and 
the Vlasov equations are conformally invariant. This implies that, provided we use conformal 
coordinates and rescaled gauge fields, the system of equations which we would have in flat 



space looks exactly the same as the one we are discussing in a curved FRW (spatially 



fiat) background [p^ ]. 

The evolution equations of the gauge fields coupled to the two Vlasov equations can be 
written as 
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V ■ E = ej d^p[f+{x,p,r]) - f_{x,p,r])], 

V X E + B' = 0, 
V-5 = 0, 

VxB-E' = J (fpv[f+ (f , p, 7]) - (f, p,7])]. (4.4) 

Now, if Sf±{x, p, rj) at the beginning of the radiation dominated epoch rj^ are known (their 
magnitude follows from our computation of the Bogoluibov coefficients), and E{x,rio) = 
B{x,riQ) = initially, the magnetic field at later times can be found from Eqs. (|4.2|) - (|4.4|) . 

This problem is solved in the Appendix B in a relaxation time approximation for the 
collision integral and under the assumption of small fluctuations in the distribution functions. 
Electric fields, as well as the charge density perturbations quickly relax to zero during the 
time given by the inverse conductivity of the plasma, a ~ nq/{{p)v), where Uq ~ is the 
density of the charged particles and {p) ~ T is their average momentum, v is the frequency of 
collisions [|. As for the magnetic fields, for the most interesting case of large scales (/c^ <^ a/r] 
the result does not depend on the frequency of collisions v and reads 

B ~ i^V X J, (4.5) 

ariq 

where a is the fine structure constant. 



V. ESTIMATES OF MAGNETIC FIELD 

We are now ready to estimate the produced current and charge density fluctuations. 
We will firstly discuss the case of magnetic fields assuming that the primordial background 



^AU quantities here are given for conformal time. For instance the curved space conductivity a 
is related to the flat space one, iTc , through a rescaUng which involves the scale factor : a = aac- 
Once we use our rescaled system n, T and a do not change for adiabatic Universe expansion, 
provided that the efl'ective number of massless degrees of freedom is constant. This is diff'erent 



from Ref. |]17| where ordinary (flat space) conductivity scales as temperature squared. 
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is of de Sitter or quasi-de Sitter type. We will then move to estimate the charge density 
fluctuations. 



A. De Sitter case 



From Eq. (|3.15|) we see that the correlation function of current density fluctuations is 



determined by the sum of the Bogoliubov coefficients. From Eq. ( p.7|) we have 



«A: + A|>p + /3p|'~ (^j^^j^sm*^y-'\kr],\-^\pr],\-\ (5.1) 



which implies 



where 



X(L) = Jd'xJ A / d'k I rf3pE^!JilMe-i^e-i^e^(^^+P^-(^^-^^) ~ 0.9 x i2nfL'. 

(5.3) 

The final result for our correlator is 

1/2 " r r(3/4)-' sJl''^ ' ^ ' 

where ~ is the typical volume of a region of size L; we take ~ r]i (larger values of 
^0 give even smaller values of the magnetic fields). 

In order to get an estimate of the magnetic field, one should specify the density of the 
charged particles and their average momentum, see ( [4.5|) . A most realistic estimate would 
be to take ~ and (p) ~ T, where T is the reheating temperature. Then we obtain for 
the gauge field fluctuations 

52 fm\-^fHi\^ 1 



In Eq. ( |5.5| ) L is the coherence scale of the magnetic field. The ratio of the magnetic field to 
the is roughly constant during the Universe expansion (if dynamo effects and the galaxy 
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collapse are discarded as well as annihilation of heavy particles) and may be taken at any 
time, provided the coherence length L is taken at the same epoch. We will take L ~ 1 Mpc 
PJT7[] at the present microwave background temperature 0, which gives LT ~ 3 x 10^^ . For 
a galactic mass perturbation (of the order of 10^^ solar masses, including dark matter) the 
typical length scale is of the order of 1.9 x (fig^^)^^^^ Mpc |2^. The gravitational collapse 
of the protogalaxies enhances the magnetic flux frozen in the primeval galactic patch by 
roughly a factor of the order of 10^. 

The obtained value for the magnetic field should be compared with the values of the 
magnetic field necessary to seed the galactic dynamo mechanism. The differential rotation 
of the galaxy introduces a parity violating term in the MHD equations (the dynamo term). 
The effect of this term exponentially amplify the seed magnetic field by a factor e'"* 0,0- In 
this amplification factor enter two numbers: the galactic age (of the order of 10 Gyrs and 
the dynamo amplification rate (F) whose estimate is rather uncertain: values of the order 
0.3-0.5 Gyr for F^^ are present in the literature Following a recent analysis [j^ we 

have that the required value for the seed field can be expressed as 



> 5 X 10"^^ for F"^ = 0.5 Gyr, 



-Bdcc 
-'dec 

> 2 X 10-2^ for F-^ = 0.3 Gyr. (5.6) 

dec 

where Tdoc is the decoupling temperature. These values have been obtained in the case of 
VLq ~ 0.3, ~ 0.7 and h = 0.65 as a fiducial set of cosmological parameters p2| , |33[| . Notice 
that in the case of a flat Universe with fig = 1 we would get values of -Bdec/^dec close to 
10-1^ 

We want now to compare these values with the parameter space described by our esti- 
mate. If we take m = 100 GeV (the lower mass limit for the Higgs boson) and if we assume 



^Equally, in agreement with the analysis in Ref. 0], one can take L ~ 100 pc at the galaxy 
formation time, which gives essentially the same value for LT. 
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Hi/Mp ^ 10 ^ (i.e. the maximal Hi compatible with microwave background anisotropies) 
we have that 

~ 5.77 X 10-^°, (5.7) 

twenty orders of magnitude smaller than the required value in order to seed the magnetic field 
of our galaxy. One can argue that by lowering the mass this estimate will improve. This is 
not the case. By taking m ~ 1 MeV (which is not realistic at all) we get B^ec/Tl^^ ~ 10~^^'^, 
still too small to be relevant. Notice, finally, that to tune Hi does not help: since it can only 
get smaller that 10~^ it can only make the value of the seed even smaller than the values we 
just mentioned. 

The most conservative estimate of the magnetic field can be obtained with the asumption 
that the number density of the charged particles in the plasma is the number of gravitation- 
ally created scalars, 

f d^k ^ 1 fHi\i^,^ (Hi 



and their average momentum is simply (p) ~ i/i ~ lOOGeV, specific for the gravitational 
production. If this is indeed the case, the value of the magnetic field B is larger by a factor 



3 



if compared with the estimate (|5.7| ). Following these considerations the obtained magnetic 
field becomes 

B . /Hi\ 1 

Even with the most optimistic numbers (i.e. m ~ 100 GeV and Hi ~ 10~®Mp ) we get 
that B /T"^ can be, at most, 10~^^ over scales relevant for the dynamo action, with a minute 
gain, with respect to the estimate of Eq. ( p. 71 ). In priciple, one can think that even if 
the produced fields are too weak to turn-on the galactic dynamo they could be of some 
relevance for other processes occurring at different epochs in the life of the Universe. For 
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instance, the electroweak horizon at the time of the electroweak phase transition [^,0 gives 



LewTe^ ~ 10^^ This would imply from Eq. (|530| ) that B^^/T^^ ~ lO'^^ At the electroweak 



epoch the smallest coherence lenght of magnetic fields is set by the diffusivity. Around the 



diffusivity scale LdifrTc„ ~ 10 Thus the obtained magnetic field can be, at most, 

i?cw ~ lO^^^ew- ■'■^ order to have sizable effects on the phase diagram of the electroweak 
transition we should have, at least, that -Bew/^ew — 0-2 HH so the produced fields are also 
too small in this context. 

B. Quasi de Sitter case 

Up to now we assumed that the primordial phase in the evolution of the Universe was 
of pure de Sitter type. In this subsection we are going to study what happens is this 
requirement is relaxed. In principle we know that during the inflationary phase the scalar 
field slightly decreases its value and, consequently, the inflationary curvature scale is not 
exactly constant but mildly decreasing as a function of time. In order to account for the 
decrease in the curvature and in the scalar field we can define the two slow-rolling parameters 

. = A = |.. (5.11) 

where denotes, in this Section, the inflaton and H = a/a is the Hubble parameter (the 
over-dot denotes differentiation with respect to the cosmic time). The important point for 
our purposes is that the slight decrease in the curvature corrects the evolution equation for 
the charged scalar. More specifically, we know that the dependence on the curvature appears 
in the mode function as a" / a. It is a simple exercise to show that 

= 2a^H\l--), (5.12) 

d ^ 

which implies that the slow-rolling corrections will make the index appearing in the Hankel 
functions slightly larger 

p = ^ + e. (5.13) 
16 



the case e = corresponds to the pure de Sitter case. Clearly, from our general expressions 
of the Bogoliubov coefficients (see Eq.( [A.16D ) an increase in p implies an increase of the 



Bogoliubov coefficients in the infra-red part of the spectrum. 

Therefore we want to estimate the magnetic fields in the case where the index v is kept 
generic. For this purpose, from Eq. ( [A. 161 ) we have that 



|«. + Al + /?,r ^ (2'^-\p - sin^ ^-)p''\k^^r'^\p^,\-'^. (5.14) 

Following the same steps as in the pure de Sitter case we can estimate that 
^ ^ n^^-^p-lYT^p) , vr N (rn^ - f ^^^^ 

T^-[ 2r4(3/4)e2 '''' sJImpJ [MpJ ^ ^ ' ^^'^^^ 

If p could be larger than 1.5 the magnetic fields would also be larger at the relevant scales. 

There are two relevant bounds on p. The first and obvious one comes from the energy 
density stored in the charged scalar modes. The energy density stored in the scalar field 
modes goes as m/c^|/5fcp. This means, in the case of generic p that the energy density (in 
critical units) is 

(m;) [jfJ (r) • 

If we take m ~ 100 GeV and p ~ 2 we see that this expression still gives a value 10~^ at 
the decoupling scale. Larger values of p could induce further anisotropies in the microwave 
background. So we will assume 1.5 < p < 2 which would be already enough to increase 
magnetic fields according to Eq. ( |5.15|) . 



The second class of bounds stems from the fact that p is connected with the slow rolling 
parameters which are constrained. In fact, the contribution to the scalar spectral index 
deduced from the COBE data could be written, in terms of the slow-rolling parameters, as 
i3l,|35| 



n = l-4e + 2A. (5.17) 

The values of e and A are different depending upon the different models of background 
evolution, namely upon the different analytical forms of the inflationary potential driving 
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inflation. This can be appreciated by writing the equations of motion of the inflaton in the 
slow-rolling approximation 

• dV 

3iJ0+ — ^0, 

Mj,H^ ~ V. (5.18) 
By using these two equations we can re-express e and A as 

Clearly the values of e and A depend upon the value of 0. For instance, we could estimate 
the value of e coinciding with the value of the field 50 e-folds before the end of inflation 
(corresponding to the moment where the large scales went out of the horizon) [^. In this 
case, for a power-law potential V ~ 0'' we have that 

^ A = ^^. (5.20) 



g + 200 g + 200 

In the case of an exponential potential of the form V = exp [60^/(pMp)] we have that 
e = X = 1/p. Consequently, from Eq. ( |5.17| ), we have that n = 1 — (2 + g)/100 for power-law 
potentials, n = 1 — 2/p for exponential potentials. The scale- invariant spectrum as it has 



been observed by the DMR experiment aboard the COBE satellite the spectral index 

can lie in the range 0.9 < n < 1.5. In order to have more magnetic fields we should increase 
p, namely we should go for large e. The variation of the spectral constrains the maximal 



value of e. So if we take 0.9 as the minimal value for n we would have from Eq. ( |5.20|) that 



the q (for a power-law potential) is g = 8. But this would imply that the maximal e is 0.03 
and our effective p will be 1.53. Too small to give relevant consequences in Eq. ( |5.15| ) for the 



magnetic fields generation. Similar conclusions could be reached in the case of power-law 
potentials. By playing with the value of p it is not possible to enhance the value obtained 
for large scale magnetic fields. 
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C. Charge density fluctuations 



On the basis of the kinetic discussion the modes which survive in the plasma are the 
transverse ones. The charge density fluctuations, being associated with the longitudinal 
modes will be dissipated quite quickly in a typical time of the order of the inverse temper- 
ature. Still it is interesting to check if the charge density fluctuations are small. From our 
expression of the Bogoliubov coefficients given in Eq. ( p^.7| ) we have that 

TT^ I J^^l I 



We can insert this last expression into Eq. ( p.l2|) and perform the integrations. The inte 



grations over x and y are trivial. The integration over the moduli over the momenta leads 
to a non trivial integral which can be exactly computed: 

roc fj ■y „ /"OO 

/ — e-" / wMwe-'"\in\i2wz = -l/2. (5.22) 
JO z Jo 

Defining then the electric charge density as 

= ^ (5.23) 



we obtain that 



n-y \MpJ 

where ~ is the photon density. This value, for a length scale corresponding to the 
horizon at decoupling, would give ng/n^ ~ 10~^^. We want to stress that the bounds 
38| , |39[| on the electric charge fiuctuations were derived by assuming that charge fiuctuations 



would induce electric fields coherent over the whole horizon. These fields would cause some 
observable effect in the microwave background so that a constraint on the charge density 
could be derived. Again, on the basis of kinetic treatment of plasma fluctuations, we can 
say that electric flelds dissipate as soon as the the plasma becomes conducting. Therefore 
the effects on the microwave background are not present. 
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VI. CONCLUSIONS 



In this paper we discussed the amphfication and the fate of the fluctuations of a charged 
scalar field in the inflationary Universe scenario. These fluctuations may lead eventually 
to the generation of some magnetic fields in the Universe. We found that the produced 
magnetic field are always much smaller than the most optimistic lower bounds required in 
order to seed the galactic dynamo mechanism. Thus, the inflationary production of charged 
scalars is unlikely to be responsible for the observed galactic magnetic fields. 
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APPENDIX A: BOGOLIUBOV COEFFICIENTS 



Defining $i and $2 as the real and imaginary part of $ as 

$ = ^($1 + ^*2), (A.l) 

we assume tfiat tlie background geometry evolves from 77 ^ —00 to 77 ^ +00 for instance, 
according to Eq. (|2.6|) . In both limits we can define a Fourier expansion of $1 and $2 in 
terms of two distinct orthonormal sets of modes. By then promoting the classical fields to 
quantum mechanical operators in the Heisenberg representation we can write, for 77 —00 

'""'^ V) = J J0j2ibpfp(v)e''^' + blj;iv)e~'n- (A.2) 



where the two sets of creation and annihilation operators (i.e. [0^,0^^] and [bp,b^^p]) are 
mutually commuting. As 77 ^ +00 $1 and $2 can be expanded in a second orthonormal set 
of modes 



V)=l J^[bp9piv)e''' + bl,g;{v)e-'n (A.3) 



Since both sets of modes are complete the old modes can be expressed in terms of the new 
ones 

fkiv) =(^k9k{v)+Pk9l{v)- (A.4) 
This transformation, once inserted back into Eq. ( |A.2|) , implies that 

ak = akak + Pla^_k- (A. 5) 

Notice that in order to preserve the scalar products in the old and new sets of orthonormal 
modes we have that the two complex numbers ak and Pk are subjected to the constraints 
lo^fcP ~ l/^fcP = 1- Exactly the same discussion applies for the field operator $2- Eq- 
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( |A.5| ) is nothing but the well known Bogoliubov-Valatin transformation; ak and (5k are the 
Bogoliubov coefficients parametrizing the mixing between positive and negative frequency 
modes. 

In order to ensure the continuity of the operators $i and $2 we have to match con- 
tinuously the old mode functions with the new ones. During the primordial phase of the 
Universe the evolution equations satisfied by the mode functions is 



^2 _ a{a + 1) ^ }J_,Vi.2a 
^2 ^2 ^ 



fk = 0, (A.6) 



where a = 1 for a pure de Sitter background fi = m/{Hiai) = mrji. Notice that with 
H we will denote the Hubble factor in cosmic time ( as usual the relation between the 
Hubble factor in conformal time, e.g. Ti and the Hubble factor in cosmic time is H = Ti/a; 
during the de Sitter epoch H ~ rj^^). Notice that if the mass of the charged scalar is not 
of Planckian magnitude, /i <^ 1. Moreover, in the limit where /i > 1 one can argue that 
the amplified fluctuations will be exponentially suppressed. Since we want to explore the 
situation where the mass of the scalar field is of electroweak order we will always be (quite 
safely) in the limit /i <^ 1. 

The exact solution of Eq. ( [A.6| ) which reduces, in the limit t] —00, to the usual 
positive frequency Minkowski space solution is given by 



/fc(^) = -^P Hl,'\-x), (A.7) 



where x = kr] and H^^'> is the first order Hankel function. Again, in the pure de Sitter case. 



p = 3/2y 1 — (4/9)/i^. Since /i is typically small, p ^ 3/2 in the pure de Sitter case. We 
denoted with p a phase factor which we choose such that 

1 ^i^i^+^P) 



p = e'i(^+'PK (A.8) 



With this choice of p we have that fk{^) ~ e ^'^'^/V2k for rj —>■ —00 



After the radiation dominated phase sets in (for rj > —rji ) the evolution equation obeyed 
by the mode functions gk{f]) is given by 
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jj'^ir] + 2r]i) 



4 

Vi 



]9k = 0. 



(A.9) 



This last equation can be easily recast in the form of a parabolic cylinder equation p6| , p7 
Defining 'j = fi/rjl we can introduce two new quantities, namely 



(A.IO) 



Consequently, Eq. ( [A.10| ) becomes 

dz^ 



+ [q+ -r]9k = 0, 



4 



(A.ll) 



which is one of the canonical forms of the parabolic cylinder equation [^,^. The exact 
solutions which reduce to positive and negative frequency modes for r] +oo are 

1 



(A.12) 



where D^r are the parabolic cylinder functions in the Whittaker's notation [27|. Notice that 
with our choice of normalizations we have, in the limit z S> \q\ and for k'^rji <C m, that 



' -e 2 ^1 



2mr] 



(A.13) 



In view of the actual calculation it is worth recalling the exact expressions of the parabolic 
cylinder functions (in the Whittaker form) in terms of confluent hypergeometric (Kummer) 
functions iFi(a,6, x): 



/ — 13. — 1. 



— i2 

U 2 



1 ,1 iq 1 iz'^ . 



•4 2'2' 2 



z{l — i) ^ / 3 3 iz^ . 

't{\ - '4)^ 2"'2' 2~ 



^4 2 / 

2 4 g 4 



1 pri^!^ 1 

r(| + i)^^^^4^ 2 '2' 2 ^ 



z(l + z) p / 3 iq 3 iz'^ 



^4 ' 2 . 



(A.14) 
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The Bogoliubov coefficients are obtained from 

fk{-Vi) = ak9k{-Vi) + Pk9*k{-Vi), 

f'ki-Vi) = akg'ki-Vi) + Pkgl'hvi), (A.15) 

which is a system of two equations in the two unknowns ak and jSk- 

By solving this system we obtain an exact expression for the Bogohubov coefficients 
which is, in general a function of two variables : fi = mrji and Xi = krji. Since /i <^ 1 we can 
expand the exact result in this limit and we obtain, in the case of a generic Bessel index p. 



A = 7re-^t{-^^|^52(xi,p)p-3 + ^^[5i(xi,p) + S^^x^.p)]^"^ + 0(pt), (A.16) 
where S\{x\^p) and S^ixiiP) contain the explicit dependence upon the Hankel's functions: 



Xx, 



If we are now specifically interested in the pure de Sitter phase we can insert the value 
p = 3/2 into Eq. ( [A.17|) . Then, we can insert the obtained expressions into Eq. ( |A.16|) and 
we obtain the wanted Bogoliubov coefficients reported in Eqs. ( |2.7| ). 

Notice that it would not be correct to use the asymptotic solutions (like the one reported 
in Eq. ( |A.13| )) in order to compute the Bogoliubov coefficients^. In fact Eq. ( |A.13| ) can be 
viewed as the a WKB-type solution of Eq. ( [A. Ill ) which can be also written as 



drf 



+ uj'Mgu = 0, uj'M = k' + m'a\ri). (A.18) 



^We disagree with the calculation of Ref. |17| where the matching has been performed using 
approximate mode functions. In our case, for small /i the leading behaviour of the Bogoliubov 
coefficient goes as . In the case of |17] it goes as an artifact of the WKB approximation. 
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By now postulating a WKB-type solution we have that 



1 








2 






kWo) _ 



By now inserting the trial solution back to Eq. (|A.18| ) we get W{r]) is specified by the 
following non-linear relation 

This equation can be solved by iteration. If we keep the lowest order we get that 

Woiv) - Mv) (A.21) 

and by using the explicit expression of the scale factor during the radiation dominated 
epoch we exactly get Eq. (|A.13|) . This solution is valid provided the corrections to the exact 
expression of W{ri) are small, namely, provided, from eq, (|A.20|) 

(A.22) 

This last inequahty, using the explicit expression of uikirf), implies that 

kY + mW { "1^^] » 1. (A.23) 
V r^i / 

Now we can see that this inequality is clearly satisfied for t] +oo. However, for rj rji, 
this inequality would imply mr]i > 1 (since krji < 1). Therefore, in order to be consistent 
with the requirement that mr]i < 1 we have to use the WKB-type solution only for large 
(positive) rj. 

APPENDIX B: VLASOV-LANDAU APPROACH TO ELECTROMAGNETIC 

FIELD FLUCTUATIONS. 

The purpose of this appendix is to give details concerning the derivation of the relation 
between the electromagnetic field fiuctuations and the initial fiuctuations in the current (or 
charge) density profile. 
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By subtracting Eqs. (|4.2| ) and ([4.3| ) we obtain the equations relating the fluctuations 
of the distributions functions of the charged particles present in the plasma to the induced 
gauge field fiuctuations: 

d d ^ dfo 

T^f{x,P,t) + v- —f{x,p,t) + 2eE-—^ = -v{p)f, 
or] ox op 

V ■ E = e j d^pf{x,p,r]), 

V X E + B' = 0, 
V-fi = 0, 

V X B - E' = f d^pvf{x,p,ri), (B.l) 



where f{x,p,ri) = 6f+{x,p,ri) — 6f-{x,p,ri) and z/(p) is a typical frequency of collisions. 

We can solve this system [^] by taking the Fourier transform of the space-dependent 
quantities and the Laplace transform of the time-dependent quantities: 



— * 


poo 

Jo J 


1 d^xe 


-''■'E{x,r]), 


= 


roo 

Jo J 


1 d^xe 


-''■'B{x,v), 


f ku) 




^ d^xe~ 





(B.2) 

We have now to specify, at the initial time, the form of the perturbed distribution function 
which can be derived from the amplification studied in the previous Section. We will call 
gj^{p) the initial profile of the distribution function. Eq. ( |B.1| ) can then be re-written as 

-g^.ip) +^{k.v- uj)fap) + 2e4^^ ■ ^ = -uf, (B.3) 

zk . E^,^ = e J f^,Jp}d'p, (B.4) 

z^-4^^ = 0, (B.5) 

4^. = -k X 4^., (B.6) 

Ml - 5)4. + -Hk-E^,J = I d'p V f^,Jp}, (B.7) 



where eq. ( [B.7D has been obtained by using Eq. ( |B.6D in the (Fourier and Laplace) trans 



formed of the last of Eqs. (|B.1|) . The Gauss constraint at r/ = implies that 
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ik-Eo{k)=e j £pg-^{p). (B.8) 

If we start, at the initial time, with a given profile of fluctuations fluctuation the Gauss 
constraint determines the initial value of the electric field. The magnetic field fluctuations 
are consistently equal to zero. 

We can now separate the electric field in its polarizations parallel and transverse to the 
direction of propagation of the fluctuation. The transverse current provides a source for the 
evolution of transverse electric field fluctuations 

^^(1 - ^)Et = e J d'phJp)v^, (B.9) 

whereas the charge fluctuations provide a source for the evolution of longitudinal electric 
field fluctuations 



ik ■ E^i- = e 



' J d'phJp). (B.IO) 
In Eqs. (|B.9|) and (|B.10|) we defined the longitudinal part of the electric field fluctuations 



and the transverse electric field as 



The solution of Eq. ( p.3| ) is given by 

i[k ■ V — uo — iv) op 
where we used that dfo/dp = vdfo/dp. The longitudinal and transverse components of the 
electric fluctuations can be obtained by inserting Eq. (|B.12|) into Eqs. (|B.9|) - (p.lCI| ) 



e 



fc-' ik e\\ J "^i^k-v-u-w) 



uj'^e_L-k^J (^k . ^ - - ijy") 



Et = I d'pv \, . ^ , (B.14) 



where ey and e±_ are, respectively, the longitudinal and transverse part of the polarization 
tensor 
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6||(A:,^) = 1 - / rf^P- ^, (B.15) 



ex(fc,^) = 1 - - / rf^p--^^^^ (B.16) 

{k-v-uj-iv) op 



Now, the general expression for the generated magnetic field is 

5, = , „\ / dMv X k]^-^ . (B.17) 

The space-time evolution of the magnetic fiuctuations can be determined by performing 
the inverse Laplace and Fourier transforms: 

B(x,T]) = I e-'^''^—^, / d^ke^'-'^lv x k] f d^p . (B.18) 

J ujh4k,uj)-k^ J V (if . k - u - tu) 

In order to perform this integral, the explicit relations for the polarization tensors should 
be given. They depend on the equilibrium distribution function fo{p), which we take to 



MP) = ^^-^ (B-19) 

where T is the equilibrium temperature, p is the modulus of the momentum and n is the 
equilibrium (thermal) density of charged particles in the plasma. The normalization is 
chosen in such a way that / d^pfolp) = Uq. 
Then we have for transverse polarization 

and for the longitudinal polarization: 

„(i..) = l + ^{2 + — 1„^^-^|. (B.21) 



^ Notice that most of our considerations can be easily extended to the case of a Bose-Einstein 
or Fermi-Dirac distribution. What is important, in our context, is the analytical structure of the 
polarization tensors and this is the same for different distributions [p^. 
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Consider now the case of very small momenta k <^uj and c<j <^ i/, relevant for long-ranged 
magnetic fields. Then the computation of the integral ( |B.18[ ) in the large time limit and 



with the use of explicit form of the transverse polarization tensor in ( |B.20D gives 0: 



T - ^ 

B{x, T]) ~ exp {-k'^r]/a)V x J, (B.22) 

where a is the plasma conductivity in the relaxation time approximation, 

2e^n„ 

a = ^ B.23 

uT ^ ' 

and initial electric current is given by 

J(f) = j d^pvg^^ip) . (B.24) 

In closing our discussion of the Vlasov equation we want to briefly comment about the 
validity of our approach. The obtained results assumed that the hnearization of the Vlasov 
equation is consistent with the physical assumptions of our problem. This is indeed the case. 
In order to safely linearize the Vlasov equation we have to make sure that the perturbed 
distribution function of the charge fluctuations is always smaller than the first order of the 
perturbative expansion (given by the distribution of Eq. p.l9| ). In other words we have to 
make sure that 

\su{x,p,v)\ « fo{p), \SfA^,p,v)\ « Up)- (b.25) 

These conditions imply that 

^4i?)«/„(p,, (B,26) 



\k ■ V — uj — iul dj 
If we now define the relativistic plasma frequency as 



u^l = (B.27) 



^For small k, the equation u!'^e±{k,uj) — k'^ = defining the poles of the inverse Laplace transform 
implies uj ~ iky a. 
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we can see that the condition expressed by Eq. ( [A.20|) can be restated, for modes k < Up, 
as \Ej^^\^ < rigT (where we essentially took the square modulus of Eq. ( [B.26| )). This 
last inequality expresses the fact that the energy density associated with the gauge field 
fluctuations should always be smaller than the critical energy density stored in radiation. 
The linear treatment of the Vlasov equation is certainly accurate provided the typical modes 
of the the field are smaller than the plasma frequency and provided the energy density in 
electric and magnetic fields is smaller than T^, i.e. the energy density stored in the radiation 
background. 
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